Inclusion of Quantum Fluctuations in Wave Packet Dynamics* 
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We discuss a method by which quantum fluctuations can be included in microscopic transport 
models based on wave packets that are not energy eigenstates. By including the next-to-leading 
order term in the cumulant expansion of the statistical weight, which corresponds to the wave 
packets having Poisson energy distributions, we obtain a much improved global description of the 
quantum statistical properties of the many-body system. In the case of atomic nuclei, exemplified 
by 12 C and 40 Ca, the standard liquid-drop results are reproduced at low temperatures and a phase 
transformation to a fragment gas occurs as the temperature is raised. The treatment can be extended 
to dynamical scenarios by means of a Langevin force emulating the transitions between the wave 
packets. The general form of the associated transport coefficients is derived and it is shown that the 
appropriate microcanonical equilibrium distribution is achieved in the course of the time evolution. 
Finally, invoking Fermi's golden rule, we derive specific expressions for the transport coefficients and 
verify that they satisfy the fluctuation-dissipation theorem. 



1. INTRODUCTION 

Microscopic simulations are employed in various fields 
of physics as tools for obtaining both of qualitative and 
quantitative insight. In these approaches, the many- 
body system is usually represented by a set of classical 
variables, even when it is basically quantal in charac- 
ter. For example, in molecular dynamics the many-body 
wave function is often represented as a (possibly antisym- 
metrized) product of parametrized single-particle wave 
packets, and the equations of motion for the parame- 
ters are then derived from a suitable variational principle. 
As a result, the quantal features of the system are not 
fully incorporated in the description. In particular, the 
system is no longer quantized and the quantum statisti- 
cal properties are not properly accounted for. In effect, 
the statistical operator exp(—f3H) is being replaced by 
the mean- field approximation to the statistical weight, 
cxp(— (3 < H >), although the wave packets are not en- 
ergy eigenstates. The wave packet parameters then be- 
have classically rather than quantally. This shortcoming 
is especially significant in scenarios where the statistical 
properties play a major role, such as complex processes, 
and, consequently, the quantitative utility of the results 
needs careful assessment. 

This issue is especially relevant in the multifragmen- 
tation processes occurring in energetic collisions between 
atomic nuclei ■ This phenomenon is of interest since 
it may provide information on the liquid-gas phase tran- 
sition of nuclear matter that is expected because the 



nuclear forces resemble those of a Van der Waals gas. 
However, interpretation of the data must rely heavily on 
model simulations, due to the complexity of the colli- 
sions P-pil . Since the constituent nucleons are fermions 
and the associated quantum statistical features persist to 
rather high temperatures, close to the transition temper- 
ature from the liquid to the gas phase, it may then be 
expected that the fragmentation pattern is influenced by 
the basic quantal nature of the nuclear system. 

The outcome of complex processes, such as the frag- 
ment multiplicity distribution in nuclear collisions, is to a 
large degree governed by the statistical properties of the 
system, even if complete equilibrium is not reached. It 
is therefore important to understand to which extent the 
quantal statistical features are accounted for in the dy- 
namical model employed. In our initial study [jTT[ , we ex- 
amined the statistical equilibrium properties of Antisym- 
metrized Molecular Dynamics (AMD) S which describes 
the time-development of the Slater-determinant of single- 
particle Gaussian wave packets, and we considered espe- 
cially non-interacting fermions in one dimension, either 
bound in a common harmonic potential or moving freely 
within an interval. Our principal conclusion was that the 
average excitation energy and the specific heat, consid- 
ered as functions of the imposed temperature, generally 
behave in a classical manner when the canonical weight 
employed is consistent with the dynamics [0,E2|. Fur- 
thermore, we found that the quantal statistical features 
could be recovered by taking account of the energy dis- 
persion of each wave packet jfl| . 
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This result may at first appear somewhat surprising, 
since the underlying wave function is antisymmetric and 
the quantum statistical properties of the nuclear system 
to a large degree arises from the fermionic nature of the 
nucleons. Indeed, that work has generated some debate 
about the statistical character of wave packet dynamics, 
such as Fermionic Molecular Dynamics (FMD) JlO| and 
AMD. In particular, Schnack and Feldmeier mW (SF) 
and Ono and Horiuchi |H (OH) have presented argu- 
ments in favor of the statistical properties of the wave 
packet systems being quantal rather than classical. One 
aim of the present paper is to elucidate this issue. We 
shall show that the energy dispersion of each wave packet 
plays a decisive role for the appearance of quantal sta- 
tistical features and that it can be included in a natural 
manner by way of the cumulant expansion of the sta- 
tistical weight. Moreover, the usual molecular dynamics 
leads to the classical statistical weight, exp(— f3 < H >), 
while the quantal weight, < exp(—(3H) >, can be well 
approximated by carrying the cumulant expansion to the 
next-to-leading order in f3. 

The modification of the statistical weight and the as- 
sociated distortion of the internal structure of the wave 
packet naturally lead to a modification of the dynam- 
ics itself. In our preceding work Jl5| , we have shown 
that the energy dispersion within each wave packet not 
only leads to the proper quantum statistical properties 
but also gives rise to a Langevin-type random force in 
the dynamics. In the field of nuclear physics, dynami- 
cal calculations with fluctuations have been applied to 
various dissipative phenomena, such as nuclear fission 
p6| and (idealized) nuclear collisions JT~T| p~9|] ■ In the first 
case, the unretained microscopic degrees of freedom pro- 
vide a heat bath that induces thermal fluctuations in the 
macroscopic degrees of freedom. In the latter case, the 
stochastic nature of the residual two-body collisions in 
the Boltzmann-Langevin model produces a statistical en- 
semble of possible dynamical histories for the one-body 
phase-space density. In contrast to those situations, the 
random force considered here has a purely quantal origin, 
as it arises from the energy dispersion of the individual 
wave packets, and it acts in addition to the statistical 
fluctuations induced by the residual scatterings. 

The quantal Langevin force arises from the compos- 
ite nature of each wave packet. This can be understood 
from the fact that an energy eigenstate, which is the ob- 
servationally relevant entity, can be expanded in terms 
of the wave packets and therefore instantaneous transi- 
tions between those wave packets are possible. Further- 
more, if the system is sufficiently complex to be ergodic, 
its time evolution will reflect the associated microcanon- 
ical ensemble. Since the wave packet has a finite energy 
dispersion, it may contribute to the phase volume even 
when its expectation value differs from the given energy. 
In the usual wave packet molecular dynamics, in which 
the equations of motion are obtained on the basis of the 



time-dependent variational principle, this kind of relax- 
ation does occur because the expectation value of the 
Hamiltonian operator is conserved and so only the states 
lying on that hypersurface are dynamically accessible. 

The principal aim of the present presentation is two- 
fold. The first is to elucidate a microscopic theory that 
provides a satisfactory description of the quantum sta- 
tistical properties of many-body systems such as nuclei, 
including the characteristic quantal behavior at low tem- 
perature and the transition from the quantum liquid to 
the fragment gas phase. We will show that this can be 
achieved by employing a simple approximate expression 
for the exact statistical weight, Wp =< exp(—f3H) >, 
which is derived by the standard cumulant expansion. 
At the same time, we illustrate the importance of the 
distortion of the wave packets caused by the canonical 
operator exp(—(3H). This distortion modifies the expec- 
tation value associated with a given wave packet and is 
the key to resolving debate between our work and those 
of SF and OH. 

The second aim is to use the insight gained from static 
scenarios to develop a corresponding dynamical theory in 
which transitions between states with different energy ex- 
pectation values are possible and which relaxes towards 
the appropriate quantal microcanonical ensemble. We 
will show that this can be achieved by introducing a suit- 
able Langevin force and employing Fermi's golden rule 
for specific expressions. The temperature can be defined 
without ambiguity and the Einstein relation between the 
drift and diffusion coefficients emerges naturally. 

This paper is organized as follows. In section || we dis- 
cuss the statistical properties of wave packet systems and 
show the importance of the energy fluctuation. In section 
^ we then apply the proposed statistical model to finite 
nuclei, specifically 12 C and 40 Ca . The quantal micro- 
canonical ensemble is then discussed in section ^ and we 
present the dynamical model with quantum fluctuations. 
Finally, section @ summarizes our studies. 



2. CANONICAL ENSEMBLE OF WAVE 
PACKETS 

In order to set the framework for the discussion, we 
start by considering a canonical ensemble of parametrized 
many-body wave packets. The key quantity governing 
the statistical properties of a quantal many-body system 
is the associated partition function, 

Z((3) ee Tr e-f> A , (1) 

where H is the many-body Hamiltonian operator and (3 is 
the inverse of the temperature of the canonical ensemble 
considered. 
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2.1. Wave packets 

A central issue is how to calculate the statistical prop- 
erties on the basis of the parametrized wave packets com- 
monly employed in microscopic transport simulations. 
While our discussion and treatment apply rather gener- 
ally, we shall work within the AMD framework, in order 
to allow concrete illustrations. As the basic wave func- 
tions we thus use Slater determinants of Gaussian wave 
packets, (r\z) ~ exp[— v(r - z/y/v) 2 } The A-body 
wave packet \Z) is then characterized by its complex 
centroid parameter vector, Z = [z\i Z2, ■ ■ ■ , Za), and 
the corresponding normalized wave packet is given by 
\Z>=\Z)/J{Z\Z). 

We assume that the wave packets provide a resolution 
of unity, as is the case in AMD. The anti-symmetrization 
operator is then represented as 



A 



dT \Z >< Z\ 



(2) 



The operator A projects onto the space of anti-symmetric 
wave functions, so it acts as unity within that space. In 
the AMD paramctrization, the measure dT is given by 



dT(Z) = dct C dZdZ = det C 



d 3A Re(Z) d 3A lm(Z) 



7T 



3.4 



(3) 



In order to carry out the 6A-dimensional integrals over 
dT, we have adopted a Metropolis sampling method, 
which guarantees that the computational effort is ex- 
pended in proportion to the relative importance of the 
individual wave packets considered. 

The action associated with a particular history Z(t) is 
given by S\Z(t)\ = J dt C(t), where the Lagrangian is 



C(t) =< Z\ih^- -H\Z > 



(4) 



The usual demand that the action be stationary then 
yields the equation of motion for the parameter vector 
Z, 



dZ i „ 
= —F 

dt h 



(•5) 



where the associated generalized force has been employed 

(6) 



dZ 



The elements of the A x A coefficient matrix C are 
spatial tensors with complex elements, 



o 2 



log (Z\Z) , 



y dzidz 



and H denotes the expectation value of the energy, 



(7) 



n =< z\h\z >-- 



(Z\H\Z) 
(Z\Z) 



(8) 



The equation of motion (||) is entirely classical and, 
if it can be brought onto canonical form, the partition 
function governing the statistical behavior of the wave 
packet parameter Z is given by the standard classical 
expression, 



J dT{Z) e - pn{Z) = J dT{Z) W$\Z) , (9) 



where W^^Z) = cxp(— (3H(Z)) is the classical statistical 
weight of a given wave packet. Such a transformation has 
been shown to exist for the Gaussian wave packets em- 
ployed in AMD , leading from the wave packet param- 
eter Z to the phase-space variable W = ^/vq-\-ip/2h^/u. 
The measure is then 



dT{Z) = dWdW = [[ 



-A- dpdq 



h 3 



(10) 



The above result @ should be contrasted with the 
exact quantal partition function (|l]) which can also be 
expressed as an integral over the wave packet parameter 
space, since the wave packets employed resolve unity (see 
eq- (§), 



Z((3) = Tr c 



-0H 



dT W (Z) 



(11) 



However, the integrand is more complicated, 

W/j(Z) =<Z\ e~ fi6 \Z > , (12) 

depending not merely on the expectation value TL =< 
Z\H\Z > but also on the energy fluctuations inherent in 
the particular wave packet, as is further discussed below. 

2.2. Statistical weight 

It is a hard task to calculate the exact statistical 
weight, Wp(Z), since it contains ^4-body operators. It 
would in fact be equivalent to solving the dynamics ex- 
actly. It is therefore useful to develop an approximate 
treatment. 

The main problem arises from the energy spread of 
the wave packets. Since a given wave packet is generally 
not an energy eigenstate, it has a distribution of energy 
eigenvalues, as is conveniently described by the associ- 
ated strength function, 

p B (Z) = < Z\S(H-E)\Z> ^ 8{H-E) . (13) 

Consequently, the variance of the energy distribution is 
generally positive, 
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a 2 E (Z) = < Z\H 2 \Z > - < Z\H\Z > 2 > . (14) 

When these energy fluctuations are sufficiently small, 
the statistical weight can be evaluated in the "mean- 
field" or "static" approximation, yielding the classical 
result, 



W| tatic (Z) 



-p<Z\fi\Zz 



-PH(Z) 



Wf{Z) 



(15) 



The statistical properties are then those of a classical sys- 
tem with the Hamiltonian function Ti{Z) =< Z\H\Z >, 
i.e. the same as those exhibited by the wave packet 
parameter Z, as discussed above. This correspondence 
remains true even for the AMD model, in which the 
fermionic nature of nucleons is incorporated by employ- 
ing antisymmetrized wave packets |l^,[llj . 

We show below how to improve the calculation of the 
statistical properties when the energy fluctuations are 
significant. For this purpose, we first note that the sta- 
tistical weight can be expressed as the norm of a wave 
function that evolves along the imaginary time direction, 



-i/W e -iM| Z > 



Wp(Z) = < Z\c- fm \Z > = < Z\e 
= <Z(|)|Z(|)>, 

where the (3 dependent wave packet satisfies 



(16) 



-H\Z{f3) > 



and, consequently, changes its norm with /3. The corre- 
sponding evolution equation for Wp is given by 



d_ 

d/3 



Wp{z) = -n (z)Wp(z) , 



with the (3 dependent energy expectation given by 

< Z(f3/2)\H\Z([3/2) > 
np{ > ~ < Z{(3/2)\Z{(3/2) > ■ 

The formal solution to eq. (|l^) is 

1-0 



Wp{Z) = exp 



dp' Hp 



(18) 



(19) 



(20) 



This expression for the statistical weight is exact, and 
when the quantal state \Z > is an energy eigenstate, 
eq. ( p0| ) is the usual statistical weight Wp = exp(—[3E), 
where E is the associated eigenvalue. The difference be- 
tween ([l5]) and ( p0| ) comes from the fact that the given 
state \Z > is not an energy eigenstate so its energy dis- 
tribution must be taken into account. 

The quantal weight (^) can be calculated once we 
know the (3 dependence of the expectation value of the 



Hamiltonian operator, Tip(Z). Generally, Tip(Z) de- 
creases as the temperature is reduced, since the distor- 
tion of the wave packet caused by the canonical factor 
exp(— j3E) is then more effective. This feature can also 
be brought out by a direct differentiation of Tip, 



< Z(f3/2)\H 2 \Z(f3/2) > 



<Z{f3/2)\Z{f3/2)> 
= -<4(Z(/3/2))<0. 



-Hi 



(21) 



Furthermore, if the given state \Z > is not orthogonal 
to the ground state, as is always true for the generalized 
coherent states usually employed, then Tip(Z) converges 
to the ground state energy E gs as (3 — > 00. 

In order to develop a suitable approximation to 
Wp(Z), we perform a first-order logarithmic expansion 
of Hp(Z), 

logHp(Z) «logtt(Z) + /3^1ogtt /3 (Z)| /3=0 , (22) 



and find 



Hp{Z) w H(Z) e - 0D ^ 



(23) 



Here the effective level spacing D is equal to the relative 
energy variance of the state considered, 



D(Z) = ~]ogHp(Z) 



*l(Z)/H(Z) . (24) 



0=o 



(17) 

Eq. (p0|) then yields the corresponding statistical weight, 



W (Z) « exp 



--(1-e-H 
D V ; 



(25) 



where the ground-state energy has been subtracted from 
the Hamiltonian, H(Z gs ) = 0. 

In order to elucidate this approximation, we consider a 
harmonic oscillator with the level spacing A and employ 
wave packets of coherent form, \Z) = exp(— Za^O >. 
The normalized state is then 



\Z >= e~i zz 



2?n 



(26) 



and its mean energy is Ti-(Z) = ZZA. Moreover, the 
corresponding spectral distribution is of Poisson form, 



H/A 



-t- (ZZ)" e- ZZ 



(27) 



The (3 evolution of the wave packet \Z > is easily ob- 
tained, 



W/2) > 



e^ ZZ 



E4f 



e -in/3A 1 > ^ 



(28) 



so, according to eq. (jig), the statistical weight then be- 
comes 
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W (Z)=e~ ZZ Y, ~ CZZYe 

n 

= exp [-ZZ(1 - e^ A )] . 



n/3A 



(29) 



The j3 dependence of TL can then be obtained by differ- 
entiation, 



H„{Z) = -^-logWp(Z) = ZZA e"^ A 



H(Z) e 



-0A 



(30) 



Thus the Hamiltonian Ti.p(Z) exhibits a pure exponen- 
tial damping. It then follows that the associated effec- 
tive level spacing, D = —d\ogH/d(3, is simply equal to 
the level spacing of the oscillator, A. Furthermore, we 
see that the expression (|2^) for the statistical weight 
in fact agrees with the exact result (|29|), with TL/D 
representing the mean number of quanta in the state, 
< n >= H/A = ZZ. 

The expressions (( |2q) a nd (Eij)) are those already used 
in our previous work . It was shown there that they 
yield the proper statistical properties, in certain exactly 
soluble cases. We also note that the statistical weight 
(Eq) can be obtained from the general cumulant expan- 



logWp(Z) = -pH(Z) + -p 2 H(Z)D(Z) 
= -m{Z)+ l -p 2 al(Z) + -- 



(31) 



The lowest-order term is recognized as the static result, 
eq. (|l5|), while the inclusion of the next term yields the 
present refined treatment, eq. (|25|). By carrying the cu- 
mulant expansion to higher order, increasingly refined 
approximations can thus be developed, if called for. The 
approximation adopted here is exact in the case of dis- 
tinguishable particles in a harmonic oscillator potential. 
Therefore we refer to this treatment as the Harmonic 
Approximation. 



2.3. Expectation values 

The effect of the canonical operator exp(— /3H) on a 
wave packet is two-fold. The most obvious effect is the 
assignment of a statistical weight for the wave packet as a 
whole, Wp(Z), as described above. However, since each 
energy component is weighted differently (namely accord- 
ing to its eigenenergy) , the statistical operator also intro- 
duces a distortion of the internal structure of the wave 
packet. This latter effect generally modifies the expecta- 
tion value of quantal operators and, as a consequence, the 
evaluation of thermal expectation values is less straight- 
forward than one might at first have thought. 

In general the canonical expectation value of a quantal 
operator O is given by 



Tr 



O e 



-0H 



(32) 



Using the resolution in terms of the wave packets, eq. (g), 
the above expectation value can be written on standard 
form as a weighted average, 



1 



dT{Z) Op{Z) Wp(Z) , (33) 



where the value of the observable associated with a given 
wave packet \Z > is the expectation value of O in the 
distorted state \Z(/3/2) >, 



O p (Z) = <O> = 



< Z\e-^/ 2 6e-f 3 "l 2 \Z > 



< Z\ c-f )H \Z > 
< Z((3/2)\6\Z(P/2) > 
< Z(0/2)\Z(fi/2) > ■ 



(34) 



When the temperature is sufficiently high in compari- 
son with the effective level spacing, T 3> D — a%/Tl, the 
distortion is very small and we have, 



ci 



f 



Z(J3) 



dV{Z) < Z\0\Z > Wp{Z) 



(35) 



It is relatively easy to exhibit the effect of the thermal 
distortion of the wave packets when the temperature is 
high. An expansion of eq. (|34|) shows that the leading 



correction to the classical result (35) is proportional to 
the quantal correlation between the observable O and the 
Hamiltonian iJ, 



<Z\ l -{6.H}\Z>-0{Z) H(Z) 



(36) 



where {•, •} denotes the anti-commutator. Thus, the ex- 
pectation value of O in the thermally distorted state de- 
creases if the correlation is positive and conversely. 

However, in cases of practical interest the thermal dis- 
tortion is not negligible and must be taken into account 
in order to obtain quantitatively useful results . The 
effect is well illustrated by the behavior of the mean en- 
ergy of the thermal ensemble of wave packets. Assuming 
that H.r(Z) exhibits a simple exponential evolution, as 
in eq. (|23|), the mean energy can be expressed as 



<n> a = 



-L J dT (Z) H(Z) e-? D Wp{Z) . (37) 



The distortion factor e~^ D suppresses the contribution 
from a given wave packet (as it should, since it is posi- 
tively correlated with itself), and it ensures that the char- 
acteristic quantal behavior, E* cx e~@ D , emerges at low 
temperature, T <C D. 
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2.4. Discussion and illustrations 



Quantal statistical weight: 



In our initial work, we pointed out that the standard 
AMD model exhibits classical statistical properties and 
suggested that the situation be remedied by including of 
the energy fluctuation of each wave packet |0]. That 
work has generated some debate about the nature of the 
statistical properties of wave packet molecular dynam- 
ics, such as FMD and AMD. In particular, Schnack and 
Feldmeier §]| (SF) and Ono and Horiuchi Q (OH) have 
made the counterclaim that the statistical properties of 
the wave packet systems are quantal rather than classi- 
cal. Hoping to elucidate the issue, we now discuss and 
illustrate the situation in some detail. 

In their analyses, both SF and OH focus on the oc- 
cupation probability for the single-particle levels in a 
harmonic-oscillator potential. Specifically, they consider 
the quantity 



-< Z\e^{-/3H)\Z > 



Pn(f3) 



1 



dT{Z) < Z\O n \Z > Wp(Z) , (38) 



where O n = <v n a n is the one-body operator counting the 
number of particles in the level n. In the work by SF, the 
integral over dT is carried out through the time-evolution 
of FMD wave functions in the harmonic oscillator poten- 
tial with a perturbative residual interaction; they have 
also studied the thermalization between two oscillators 
with different frequencies. The work by OH employs the 
same wave packets and measure as our original work [ p"T| , 
and they have also studied actual nuclear systems by us- 
ing the evaporated nucleons as a thermometer. The re- 
sults of these studies of the harmonic oscillator can be 
summarized as follows. When the initial energy (SF) or 
the temperature (OH) is chosen so as to reproduce the 
mean energy of the corresponding quantal canonical en- 
semble, then the quantity P n emerges as being very simi- 
lar to the occupation probability of the quantal canonical 
ensemble. 

This result is not unexpected, since both of FMD and 
AMD are based on totally anti-symmetrized wave func- 
tions and therefore take account of the Fermi statistics 
governing the individual particles. Moreover, the rela- 
tion between temperature and energy in the ensembles 
employed by SF and OH is by construction the quan- 
tal one. However, this type of analysis does not address 
the key issue, namely the many-body properties of the 
system. In order to clarify the situation, we discuss the 
following three separate issues: 

1. Statistical weight. The statistical weight Wp em- 
ployed in the integration over parameter space can 
be calculated either quantally (Q) or classically 
(CI), 



• Classical statistical weight: 
Wf = exp (-/3H) 



(39) 



(40) 



2. Expectation value. The value of the observable as- 
sociated with a given wave packet can be calculated 
with either the thermally distorted states (D) or the 
undistorted states (U), 



• Observation with distorted states: 



< Z|e-^/ 2 Oe-^/ 2 |Z > 
<Z\e-P 6 \Z> 



(41) 



• Observation with undistorted states: 



C¥ =< Z\6\Z > = 0{Z) (42) 



3. Temperature. The temperature can be either taken 
as the specified one (T) or readjusted so the expec- 
tation value of the energy is matched exactly (T 1 ), 

• T: The given temperature T — is the one 
used for the statistical weight and the obser- 
vation. 

• T'\ The given temperature T is regarded as 
spurious and is replaced by T" which is ob- 
tained by demanding that the mean energy of 
the system match the exact value. 

The exact calculation discussed in the previous section 
then corresponds to the option [Q-D-T] , whereas the cal- 
culation by OH corresponds to [Cl-U-T']. Although SF 
perform the ensemble sampling by means of the time evo- 
lution, their residual interaction ensures ergodicity and 
the treatment by SF is then be essentially the same as 
that by OH. 

To illustrate the effect of the different treatments, we 
consider four fermions in a harmonic oscillator and cal- 
culate the occupation probability of each single particle 
level, for various temperatures. Figure [j] shows the re- 
sults based on either distorted (D) or undistorted (U) 
wave packets and employing the options [Q-T], [Cl-T], 
and [Cl-T'] in each case. The quantal statistical weights 
were calculated in the harmonic approximation, and the 
quantal observation is carried out by assuming the (3 evo- 
lution of Zi to be Zj(/3/2) = Zj(0) exp(— fiTiw/2). 
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Distorted States [D] 



Un-Distorted States [U] 




FIG. 1. Occupation probability in a harmonic oscillator. 

The occupation probability for four identical fermions in a harmonic potential calculated with either distorted 
states (D: left) or undistorted states (U: right) for the observation of the occupation number and employing the 
options [Q-T], [Cl-T], and [Cl-T'] in each case. About 10 4 states were sampled for each temperature by means of 
the Metropolis method. The exact quantal results are shown by the solid dots. 



Comparison with the exact result (solid dots) shows 
that both [Q-D-T] and [Cl-U-T'] reproduce the exact re- 
sults very well, which is in accordance with the reports 
by SF and OH. Even though the latter treatment is not 
formally correct, this agreement is not unexpected, as 
we shall now discuss. The key to achieving a correspon- 
dence between the two treatments is to regard the dis- 
torted state |Z(/3/2) > as the state \Z' > sampled by 
SF and OH. Then, in the harmonic approximation where 
Hp = Hcxp(—/3D), the quantal statistical weight VVg 
can be rewritten in terms of a modified inverse tempera- 



ture P' and the expectation value with respect to \Z' >, 
logW?(Z) S -^(l-e-<">) 



nzy° {i _ e _, D) 



-0"H(Z') = \ogW2(Z') , (43) 



where the modified temperature is determined by 

e tm -1 1 
? = — =P(1 + -PD + ...). 



(44) 
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Consequently, the quantal partition function can be ex- 
pressed in terms of the classical statistical weight for the 
distorted state, 



Z Q (f3) = J dT{Z) W®(Z) = J J dT(Z') W^(Z') 



(45) 



where J is the Jacobian associated with the transforma- 
tion from Z to Z' , 



J 



w 


w 


detC 


Z 


z 


w 


w 


detC' 


z' 


z' 



(46) 



Comparing the above reformulation (|45| ) with the classi- 
cal partition function, 



z c \(3') = J ar(z') w^(z') 



(47) 



we see that if the Jacobian J is independent of Z then 
the two partition functions are proportional, 



Z^(f3) = J Z c \[3>) , 



(48) 



In addition, the expectation value Op(Z) is then the 
same as < Z'\(J\Z' > by definition. Consequently, the 
thermal averages of the observable O are also equal, 



■< U ^ " Z Q {/3) 



dv{z) o9(z) w?(z) 



j J dv(z') o°> w$(z') 



J Z Gl (f3) 



(49) 



This result brings out very clearly how it is possible, for 
each given temperature T, to introduce a modified tem- 
perature T' so that the quantal result is reproduced. 

However, this kind of treatment leads to inconsistent 
thermodynamics. The mean energy -< TL >~ should be 
given by —dlogZ/d/3, and that is indeed the case in the 
approximate quantal treatment of ref. [ fi"H . But the mod- 
ified treatment in [Cl-U-T"] yields the mean energy 



-< n yf = 



^ I dT(Z') H(Z') W°}(Z>) , (50) 



Z cl {/3') 



while a differentiation of the partition function ( |45| ) leads 
to 



dp 



— l»g^ cl = J dT(Z') H(Z') W$}(Z') , 



(51) 



which contains the additional factor d(3' /dfi = exp((3D). 
Therefore, the modification of the temperature parame- 
ter renders the thermodynamics inconsistent. 



A further difference between the two treatments is the 
partition function itself. The classical partition func- 
tion Z cl deviates from the quantal partition function 
Z® by the Jacobian factor J. If the determinant ratio 
detC/ detC is ignored, we have J ~ exp(— Aj3D). This 
ideal situation is realized in the case of A distinguishable 
particles in a harmonic oscillator potential, where the 
level spacing D = Tito is a constant. Furthermore, the 
Jacobian J is independent of Z but still depends on /?, 
of course. The partition functions Z^ and Z C1 are shown 
in fig. ^|. While the quantal partition function approaches 
unity when T — > 0, the classical partition function tends 
to zero. This occurs because the phase space for Z' is 
drastically reduced relative to the original phase space of 
Z when the temperature is smaller than the level spac- 
ing, as already discussed. The fact that the phase space, 
and hence also the statistical fluctuations, are too small 
in the ensembles employed by SF and OH may present 
a practically important shortcoming of those treatments, 
since the partition function is the key quantity governing 
the fragment size distribution in actual nuclear collision 
processes. 




0.5 1 1.5 2 
Temperature 



FIG. 2. Partition function for a harmonic oscillator. 
The partition function for A distinguishable particles 
in a harmonic potential, calculated either quantally 
(Q: solid curve) or classically (CI: dashed curve). The 
temperature is expressed in units of the oscillator spac- 
ing D — Tiuj. 



3. STATISTICAL PROPERTIES OF FINITE 
NUCLEI 

In the previous section, we employed a few simple sol- 
uble cases to illustrate the general treatment. We now 
turn to more realistic scenarios and consider the statisti- 
cal properties of finite nuclei. 

As we have seen, we need to determine the (3 evolution 
of the wave packet. In the harmonic approximation, this 
is equivalent to calculating the energy dispersion u\(Z) 
of the wave packet, since level spacing D = <7 2 E /TL is the 
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key quantity determining the statistical weight. More 
generally, the equations of motion (EOM) for the wave 
packet parameters Z can be employed for this task by 
replacing the real time t by the imaginary time —ih(3. In 
the AMD model, 



E and C v of 12 C 



dZ_ 
~dt 



-ih/3 



dZ_ 



(52) 



The first equation governs the usual time development, 
while the second one can be used for the estimation of 
the statistical weight. It may be noted that the second 
equation is usually employed for the construction of the 
ground-state nuclei, such as those used as initial con- 
ditions in nucleus-nucleus collisions, since it effectively 
cools the system down. It is therefore refered to as the 
"cooling" equation |s|. 

In the present work, we have adopted eq. ( p5| ) as the 
statistical weight, 



W P (Z) = exp 



D y ' 



(53) 



and the effective level spacing associated with a given 
wave packet | Z > is calculated by the AMD cooling equa- 
tion, 



1 

m 



/3=0 

dU dZ 
dZ~d/3 4 



dU dZ 
dZ~df3 



- C 
= H 



(54) 



Working in the harmonic approximation where = 
Hexp(—f3D), we wish to calculate the temperature de- 
pendence of the mean energy -< Ti yp and the specific 
heat, 



Cv(P) 



= (3 2 U < H 2 > y 



<H>-' 



(55) 



where the expectation value of H 2 is given by 



<H 2 >!3= D U{Z) e-? D = a 2 E {Z) . (56) 

The harmonic approximation demands that D is constant 
along the (3 evolution described by the cooling equation. 
This is not strictly true for real nuclei, since the energy 
surface is not exactly quadratic. However, at low excita- 
tions a quadratic behavior holds approximately since the 
first derivative of TL vanishes in the ground state. More- 
over, at high temperatures where (3 expansion is valid, 
the quantum fluctuaton effects are well described by the 
cumulant expansion and the harmonic approximation is 
then reliable in this regime as well. Therefore we ex- 
pect the harmonic approximation to provide a reasonable 
starting point. 








AMD 

LD with Finite T 
Free Gas 



5 10 15 20 

Temperature (MeV) 

FIG. 3. Excitation energy and specific heat of 12 C. 



The mean excitation energy per nucleon E* /A divided 
by the temperature T (top) and the specific heat per 
particle Cv/A (bottom) for a canonical ensemble of 
12 C nuclei calculated with the AMD model, as a 
function of the temperature T. Squares, triangles and 
circles show the results obtained with various values 
of the radius constant, ro = 1.2,1.5,2.0 fm, respec- 
tively. Also shown are the results for a free classical 
gas (dashed lines) and a nuclear liquid drop model 
(solid curves) in which the excitation spectrum con- 
tains the known low-energy levels plus a Fermi-Dirac 
gas of quasi-particles with a suitably modified level- 
density parameter, a = A/{8 MeV)(l - 0.8/A 1/3 ) @. 
The marks represent results with three different freeze- 
out volumes, V = 4nR 3 /3,R = r A 1/3 , r =1.2, 1.5, 
2.0 fm. The sample size for each scenario is 5 x 10 4 . 
In order to give an impression of the sampling error, 
the results of two different ensembles are shown for 
r = 2.0 fm. 



Using the above framework, we have studied the statis- 
tical properties of 12 C and 40 Ca with the AMD model. 
The effective nuclear interaction used is Volkov No. 1 pCfl 
and the width parameter v and the Majorana mixture are 
chosen so as to fit the binding energy and the r.m.s. ra- 
dius. The phase space integral are evaluated by means 
of a Metropolis sampling method. 

Figure || shows the excitation energy E* and the spe- 
cific heat Cv of 12 C as a function of the temperature. 
For convenience, the excitation energy has been divided 
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by the temperature (the free value is then a constant 
(= 3/2)). For comparison, we also show E* and Cy 
for free classical particles (dashed line) and for a finite- 
temperature liquid drop model that includes known low- 
energy nuclear levels augmented by a Fcrmi-Dirac gas 
(solid curve). 

The calculated results exhibit the expected behav- 
ior: At low temperature the excitation energy behaves 
quadratically, E* j A rs aT 2 , and as the temperature in- 
creases it approaches its free value, Ef lee — 3T/2. This 
behavior is consistent with the transition from the liq- 
uid phase (or quantal phase) to the gas phase (or clas- 
sical phase). The same features appear for the specific 
heat: It grows approximately linearly at low tempera- 
tures, Cy/A w 2aT, and then approaches its free value, 
Cy ee /A = 3/2. In addition, the low-temperature results 
obtained for normal density, ro = 1.2 fm, are very close 
to those of the standard liquid drop model. These results 
are very satisfactory, since they yield a good description 
of the statistical properties of finite nucleus over the en- 
tire range of temperatures. 

Specific Heat of 12 C 




Temperature (MeV) 

FIG. 4. Specific heat and fragmentation of 12 C. 

The specific heat per nucleon calculated with AMD 
(solid dots) and with the canonical multifragmenta- 
tion model described in the text admitting a maximum 
of N fragments, with N=2-9 (short dashes). The ra- 
dius constant was taken as ro = 2.0 fm, corresponding 
to a freeze-out density equal to about 20% of normal. 
Also shown are the results for the liquid drop model 
(solid curve) and the free nucleon gas (long dashes). 

It is interesting to study the volume dependence of 
these quantities, since the freeze-out volume is the most 
important parameter in the statistical models of multi- 
fragmentation, i.e., it determines the ratio between bi- 
nary and multifragment decays pl[. The qualitative be- 
havior (the transition from the liquid-like phase to the 
gas-like phase) does not change with volume, but it be- 
comes easier to excite the system when the volume grows 
larger. This is because the excitation of intrinsic modes 



in a single nucleus is being overwhelmed by the agita- 
tion of multi-fragment configurations. In addition, there 
appears a region in which the specific heat exceeds its 
free value of 3/2, reflecting the rapid activation of many 
degrees of freedom that were effectively frozen at low 
temperatures. Thus, such a rapid increase of the specific 
heat may signal the onset of multifragmentation. 

In order to explore the validity of this picture, we have 
performed a simple statistical calculation on the fragment 
yields. The system is assumed to appear as a configura- 
tion of distinct exitable nuclear fragments, {rt;}, where 
rii is the multiplicity of the nuclear species i. The total 
partition function is then given by 

2 = }EII^ Z T G ~ mM) . ( 57 ) 

where Z t is the contribution to partition function arising 
from a fragment of type i within the volume 17, 

*- n (ssO c " ,58) 

and Z refers to the compound nucleus Z. The intrin- 
sic degrees of freedom are taken into account through the 
effective intrinsic partition function Q which contains the 
ground-state degeneracy gi — 2Jf' s ' + 1 and the excited 
levels up to E* — Eq + A (Ai — 4) for fragments with 
Aj > 4, where Eq is the minimum energy for particle 
evaporation and A = 2 MeV, as suggested in ref. p2[ . 
Furthermore, the potential energy includes the binding 
energy difference and the inter-fragment Coulomb poten- 
tial, 

V{{n i }) = B{A,Z)-Y J B{A a ,Z a )+Y j y a fi , (59) 

a a</3 

where a and (3 enumerate the fragments in the particu- 
lar configuration {r^}. The fragment binding energies B 
are taken from experiment and the Coulomb interaction 
is fixed by the barrier height. Figure |^ shows the re- 
sulting the specific heat per nucleon calculated with this 
fragment canonical model with up to nine-fragment con- 
figurations included, as well as the AMD results from fig. 
^ with r$ — 2.0 fm. It is seen that the AMD behavior is 
well reproduced when a sufficient fragment multiplicity 
is admitted in the statistical calculation. It is especially 
noteworthy that the reproduction of the rapid increase at 
T ~ 4 MeV, requires the inclusion of channels with four 
or five fragments. This illustration shows that multifrag- 
mentation appears naturally as the volume is increased 
and, moreover, the model presented in this work can de- 
scribe this phenomenon. 

It is necessary to study statistical properties of heav- 
ier nuclei in order to fully elucidate the relation between 
multifragmentation and the liquid-gas phase transition, 
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since the latter is closely related to the properties of nu- 
clear matter. Figure || shows the calculated excitation 
energy and specific heat for 40 Ca . Since the sampling 
number at each temperature and volume is only around 
2 x 10 3 , the statistical error is not small. However, the 
qualitative features are almost the same as for 12 C and, 
in particular, the low-temperature correspondence with 
the finite-temperature nuclear liquid-drop model is good. 
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FIG. 5. Excitation energy and specific heat of 40 Ca. 



The mean excitation energy per nucleon E* /A divided 
by the temperature T (top) and the specific heat per 
particle Cv/A (bottom) for a canonical ensemble of 
40 Ca nuclei calculated with the AMD model, in a dis- 
play similar to fig. ^. The level-density parameter is 
A/(8 MeV). 



DYNAMICS WITH QUANTAL 
FLUCTUATIONS 



Until now we have studied the statistical properties of 
an assembly of nucleons. Below we briefly describe how 
the insights gained can be applied in dynamical scenar- 
ios. The most important result of the previous analyzes is 
that the quantal nature of the system demands that the 
energy dispersion of each wave packet be taken into ac- 
count. This follows from the fact that the wave packet is 
not an eigenstate of the Hamiltonian operator and there- 
fore the observed energy exhibits fluctuations around the 
given initial energy value jL5|. This feature then affects 



the dynamical evolution significantly, as can be most eas- 
ily seen from the expression for the microcanonical phase 
volume, as shown below. Since the long-time dynamical 
evolution of a closed system should reflect the associated 
microcanonical ensemble, we first discuss the statistical 
properties of a microcanonical ensemble. 



4.1. Quantal microcanonical ensemble 

For an ergodic system, in which each state is reach- 
able from any given initial state, the long-time dynam- 
ical evolution should reflect the associated microcanon- 
ical ensemble. Therefore, the statistical properties of a 
dynamical model are related more closely to the micro- 
canonical ensemble than to the canonical ensemble. The 
microcanonical phase volume is given by 



Q(E) = Tr 5(H - E) 



dT 



< Z\6{H~E)\Z > 



< Z\Z > 



J dT p E (Z) 



(60) 



Since the wave packets are not energy eigenstates, i.e. 
p E {Z) 7^ 5(H — E), those states with H ^ E can also 
contribute to the microcanonical phase volume accord- 
ing to the probability of p E (Z). It is interesting to note 
that resolution of unity (|2|) causes the strength function 
p E {Z) plays a dual role: While it is defined as the energy 
distribution in a single wave packet Z ', it also expressed 
the relative weight of different wave packets in the mi- 
crocanonical ensemble characterized by the energy E. 

From the microcanonical phase volume (^0|), the en- 
semble temperature T can be defined unambiguously by 
means of the formal thermodynamical relation as the in- 
verse of the rate at which the phase volume increases 
with energy, 



hjE lo ^ = m 



dp E (Z) 
dE 



1 



fl(E) 



J drp E (z) Pe (z), 



(61) 



where the contribution from a given wave packet is given 
by 



(3 B {Z) = A log Pe (Z) 



(62) 



We note that this latter state-dependent quantity can 
have either sign. (Since p E (Z) is largest when E w H, it 
will generally increase with E when E < TL and decrease 
when E > H.) 

The microcanonical phase volume (|6^) is intimately re- 
lated to the partition function of the associated canonical 
ensemble, 
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-PE 



rOO p pQO 

Z((3) = dE Q(E) e~ pE = dT dE p E {Z)e 
Jo J Jo 

dT Wp(Z) . (63) 



Thus statistical properties of both the microcanonical 
and the canonical ensemble can be studied on the same 
footing through p s (Z) and Wp(Z). We also note that 
when the energy eigenvalue distribution of the given state 
\Z > is of Poisson form, we recover the same statistical 
weight Wp(Z) as found earlier (eq. (p3|)). 
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FIG. 6. Ensemble sampling of microcanonical temperature. 

The inverse temperature of distinguishable particles 
in a harmonic oscillator as a function of the energy 
per particle, as obtained by performing a Metropolis 
sampling of the state dependent temperature ([ii]) us- 
ing the probability (^) for N =4 (squares) and 10 
(diamonds) (10 4 states were sampled for each given 
energy). Also shown are the exact analytical results 
for the corresponding microcanonical (solid), canoni- 
cal (dashed), and classical (dotted) ensembles. 

In order to illustrate the situation, we consider a mi- 
crocanonical ensemble of N distinguishable particles in 
a common harmonic oscillator and calculate the corre- 
sponding inverse ensemble temperature by means of eq. 
(|6l|). The energy distribution p E (Z) is assumed to be a 
continuous Poisson distribution, 



p E {Z) oc 



El 



T{E + l) ' 



(64) 



using units in which Tim = 1. (The number of states for a 
given energy is then Q(E) = T(E + N)/[T(E +1)T(N)}.) 
Using the above expression as the weight function, it is 
then possible to employ the Metropolis sampling tech- 
nique to calculate the average in eq. (|6l|), -< (3e >■ The 
result is shown in fig. [6] and leads to a perfect reproduc- 
tion of the exact inverse temperature for the microcanon- 
ical ensemble, as one should expect. The temperature 
depends on the number of particles N, even though the 



particles are distinguishable, and may be compared with 
the results for the corresponding canonical and classical 
ensembles. 



4.2. Quantal Langevin force 

It follows from the previous discussion that after the 
system has undergone a sufficiently long time evolution 
the probability distribution <fi(Z) for finding the state 
\Z > should be proportional to the corresponding level 
density, 



P E (Z) 



= c -nZ) 



where 



T=-]agp B (Z) 



(65) 



(66) 



plays the role of an effective potential for the param- 
eter Z. The evolution of the wave packet parameter 
vector, Z(t), as determined by the equation of motion 
(H), is entirely deterministic, without any physical ef- 
fect of the spectral structure of the wave packet. Conse- 
quently, as we have discussed above, the system does not 
relax towards the appropriate quantum statistical equi- 
librium. Furthermore, this malady is inherent in the clas- 
sical treatment and it could not be remedied by merely 
including a collision term in the standard manner, since 
the evolving wave packet would then still remain on the 
cqui-Hamiltonian surface TL — E and, therefore, the dis- 
tribution 4>{Z) would not relax to the microcanonical en- 
semble (p^). 

In order to provide the system with an opportunity 
for exploring and exploiting the various eigencomponents 
contributing to its wave packet, we wish to augment the 
equation of motion (||) by a stochastic term that may 
cause occasional transitions between different wave pack- 
ets. The analogy of the distribution ( |65| ) with the clas- 
sical Boltzmann distribution Wp(Z) = exp(—(3H) sug- 
gests a way to introduce such a stochastic term, since 
the latter distribution can be obtained by introducing a 
Langevin term in the equation of motion. Analogously, 
the quantal distribution can be obtained by augmenting 
the deterministic equation of motion (||) by a Langevin 
term & 



h 



g-C, 



(67) 



where the dot over Z indicates the rate of change over 
and above that given by eq. (||) and the random complex 
vector £ represents white noise, 



^„WO(0^ 2<W 6(t-t'). 
In order to see that the Langevin equation 



(68) 
indeed 



can produce the desired equilibrium distribution (65), it 
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is helpful to consider the corresponding Fokker-Planck 
equation, 



Zi,- ■ ■ , z A ; t) 

' A a 



d 2 



dz n dz n > 



M nn , + c.c. 



(69) 



The transport coefficients V and M govern the early 
growth rate of the average value of the parameter vector 
and the associated covariance tensor, respectively, start- 
ing from a given value Z. Therefore, the diffusion coeffi- 
cient is given by M = g ■ g. Furthermore, if we demand 
that the quantal distribution (|6^), <fi ~ cx p( — •^ r )i be an 
equilibrium solution to this equation, the form of the drift 
coefficient follows, 



dT 

dz n , 



E 



d 



dz n . 



(70) 



The second term represents the noise- induced drift which 
enters in the general case when the tensor M depends on 
Z p3| |. Since this dependence is often unimportant, we 
shall ignore this term in the following, for convenience. 

Since T is predominantly determined by 7i, its deriva- 
tive is proportional to the driving force dH/dZ and the 
drift coefficient can then be rewritten, 



" < 71 > 



We have here introduced the inverse state-dependent 
canonical temperature, 



(72) 



and the relation (|71j) is recognized as a manifestation of 
the Einstein relation. 

The two state-dependent temperatures introduced, 
Pe = —dT/dE and (3-n = dT '/OH, are identical if p E (Z) 
depends on E and Z only through either H — E or H/E. 
This is approximately the case in most cases of interest 
and they can therefore be regarded as being very similar. 
For exa mple, for the distinguishable particles considered 
in sect. 4 4.1 , with the strength function given by (64), 
we have 



0H 



n-E 
n 



(73) 



/3 E =log(W)- — lo g r(£+l) 

= hi + + • • • . 



iog(n /e) 



(74) 
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FIG. 7. Temporal sampling of microcanonical temperature. 

The inverse temperature of distinguishable particles in 
a harmonic oscillator as a function of the energy per 
particle, as obtained by performing a sampling over 
the stochastic time evolution of a single initial state. 
The systems considered and the display are the same 
as in fig. |^ and 10 4 states were sampled for each given 
energy as well. 



In order to illustrate the utility of the stochastic treat- 
ment, we show in fig. [j] the inverse ensemble temperature 
as obtained with the proposed Langevin equation (J6J])- 
The systems are the same as those considered in fig. |6|and 
the integral over states is calculated by averaging over the 
configurations obtained at 10 4 time steps, starting from 
a single random initial wave packet. The results demon- 
strate that the quantal microcanonical distribution can 
be realized as a result of the stochastic pseudodynamics 
produced by eq. (|6^) . This fact can be also seen from the 
resulting distribution of energy expectation values P(Tl), 
as is illustrated in fig. || for N=10. Since the probability 
for finding the system with a given value of Z is propor- 
tional to p E (Z), the Hamiltonian distribution should be 
given by 



P{H) dH 



PE (z) ^ m 



Pe (Z) H*- 1 dH . (75) 



The results indicate that this is indeed the case. 

Although the stochastic term appearing here is some- 
what similar to the usual random force that arises from 
the interaction between the particles under consideration 
and the heat bath in which they are embedded, it is im- 
portant to recognize that the present random force arises 
from the energy dispersion of each wave packet and thus 
has a purely quantal origin. Therefore, we refer to this 
random force as the quantal Langevin force. We show be- 
low how this quantal Langevin force can be implemented 
once the interaction between the particles has been spec- 
ified. 
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FIG. 8. Microcanonical distribution of the Hamiltonian. 

The distribution of the expectation value of Hamilto- 
nian operator, H, for the systems considered in fig. Q 
for the case with N=10 particles, based on samples of 
either 10 4 (squares) or 3 x 10 5 (diamonds) consequtive 
states generated by the proposed stochastic dynamics. 
Also shown are the result for the corresponding micro- 
canonical ensembles (dashed curve) and that obtained 
without the phase-space factor Ti^ -1 (dotted curve). 



4.3. Characterization of the quantal Langevin force 

In the preceding, we have introduced the Langevin 
equation and shown that it leads to the desired quantal 
equilibrium distribution. However, in order to apply the 
treatment to actual dynamical processes, it is necessary 
to determine the specific properties of the transport coef- 
ficients. For this purpose, we seek guidance from Fermi's 
golden rule and assume that the differential rate of tran- 
sitions from a given wave packet \Z > to others near 
\Z' > is of the following form, 

w(Z - Z') = — |< Z' E \V\Z E > | 2 p E {Z') . (76) 

Here the operator V represents a suitable "residual" in- 
teraction and E is a specified energy which is usually 
taken as the expectation value of the originally specified 
initial state. The utility of the above form (76) is un- 
derscored by the fact that it leads to the correct micro- 
canonical equilibrium distribution. This can be verified 
by considering the equilibrium condition, 



(Z) w{Z -> Z') = 4>{Z') w(Z' -> Z) 



(77) 



which is evidently satisfied when <fi(Z) ~ pe(Z), due to 
the symmetry of the matrix element in (f76|). 

In analogy with the distortion occurring for the wave 
packets in a canonical ensemble (see sect. ||)), the initial 
and final states appearing in the above matrix element, 
\Z R > and \Z' >, are related to the actual states, \Z > 



and \Z' >, by an appropriate distortion arising from the 
microcanonical constraint on the dynamical system. In 
principle, they should be the respective components of 
the wave packets associated with the energy eigenvalue 
E. However, since the projection to energy eigenstates 
poses a hard task, we have adopted the canonical distor- 
tion as a suitable approximation, 



\Z E > 



\Z{Pn/2) > 



e -/3„g/2| Z> 

< Z\e-PnH\z > 



(78) 



Application of the cooling equation (£>2|) then yields the 
parameter characterizing the distorted wave packet, 



Z E = Z 



1 



-Jn f 



(79) 



It is easily seen that the energy expectation value of the 
distorted wave packet \Z e > is equal to E, through first 
order in (3%. 

Simple approximate expressions for the moments of 
w(Z — ► Z') can be derived by expanding the transi- 
tion rate (|7^) around Z E , where both the matrix element 
and the energy spectrum have the largest values (see Ap- 
pendix [X|) . Then the transition rate from Z to a specific 
final state Z' is given by 



w(Z -» Z') 



2tt 



|V| 2 R E {Z E ) 



x e 



-(6Z-%PnF)-C-(&Z-\f}nF) 



(80) 



where SZ = Z' - Z and V =< Z E \V\Z E > denotes the 
expectation value of the interaction in the initial state Z. 
The total rate of transitions from the given state Z into 
any final state Z' can now be readily calculated, 



/9tt 
dT' w(Z^Z') « — \V\ 2 p E (Z E ) 



(81) 



The expected number of transitions taking place during 
a small time interval At is then no = woAt, which may 
also be interpreted as the probability that a transition 
occurs during At. Therefore this total rate Wg can be 
regarded as the inverse lifetime of the wave packet \Z >. 

Since the transport coefficients V and M govern the 
early growth rate of the average value of the parameter 
vector Z and the associated covariance tensor, respec- 
tively, they can be obtained from the moments of the 
microscopic transition rate, 



V n = J dT' Sz n w(Z -> Z') 
_! dU 



-Wofl-H C 



dZ 



« w (C ) nn > , 



(82) 

(83) 
(84) 
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to leading order in (3u- We note that the Einstein re- 
lation J7l| ) is indeed satisfied by these expressions and, 
moreover, the drift coefficient is proportional to the gen- 
eralized force F, as would be expected. It may also be 
remarked that both the center-of-mass position and the 
total momentum remain unchanged on the average, since 
Xln Vn = 0> while the individual dynamical evolutions 
will exhibit diffusive Brownian-type excursions around 
those averages. This behavior is to be expected, since 
the the energy H. is no longer a constant of motion but 
will fluctuate around the specified value E. 

The above simple approximate expressions (pi]), (82), 



and (84) make it possible to simulate the Langevin evolu- 



tion by picking the stochastic changes 5z n in accordance 
with the transport coefficients V and M. It is here im- 
portant that the overlap matrix C is positive definite so 
its square root exists. Thus it is fairly straightforward to 
implement the proposed stochastic treatment. 



5. SUMMARY 

In the present work, we have addressed the treat- 
ment of quantum fluctuations in microscopic descriptions 
based on wave packets, a problem encountered in a broad 
range of fields involving quantum physics. Since the wave 
packets are not energy eigenstates, the statistical opera- 
tor exp(— (3H ) cannot be treated as a c-number. In order 
to take account of the associated spectral distribution 
and ensure that the statistical properties are quantal, it is 
necessary to introduce suitable modifications relative to 
the ordinary equations of motion for the wave packet pa- 
rameters which are basically classical in character, having 
been derived from the time-dependent variational princi- 
ple. 

We have formulated a simple but apparently success- 
ful treatment by including the first correction term in the 
cumulant expansion of the statistical weight. The associ- 
ated small parameter is a^/TH, where TC =< Z\H\Z > 
is the mean energy of a wave packet and o\ is the corre- 
sponding energy variance. This approach is exact when 
the spectral distribution is of Poisson form, as is often the 
case (at least approximately), and the corresponding ef- 
fective level spacing is D = — dlogl-tp /d [3 = a%/TL. It is 
then straightforward to write down the improved expres- 
sion for the statistical weight. Moreover, the associated 
thermal distortion of the internal structure of the wave 
packet was determined. 

Since our initial suggestion |ll|] that this treatment 
might be useful has led to some debate we dis- 

cussed and illustrated the various possible approaches to 
determining the statistical behavior of one-body observ- 
ables, such as the occupation number. In this manner, 
it was brought out that although it is possible to recover 
the quantal appearance of single-particle observables by 
suitable redefinition of the temperature parameter, such 



attempts do not yield the proper many-body properties, 
as governed by the behavior of the partition function. 
The key to resolving the issue lies in the inevitable distor- 
tion of the wave packet caused by the canonical operator 
exp(-/9ff). 

To illustrate the practical utility of the treatment, we 
considered the statistical properties of finite nuclei, as ex- 
emplified by 12 C and 40 Ca. The low-temperature behav- 
ior matches well with the finite-temperature liquid-drop 
model. Moreover, as the temperature is raised, that the 
nuclear liquid drop evolves into a fragment gas, as is ex- 
pected for actual nuclar systems and in good agreement 
with the standard statistical multifragmentation model. 

We then turned to the important issue of how to in- 
corporate the quantum fluctuations into the dynamics. 
Here the key suggestion is to allow the system to explore 
its spectral distribution by introducing suitable stochas- 
tic transitions between the wave packets fl5[| . This can 
conveniently be done by means of a Langevin term in 
the equations of motion for the wave packet parameters. 
We derived the general form of the associated transport 
coefficients and verified that the proper microcanonical 
equilibrium distribution is indeed achieved. Simple ap- 
proximate expressions for the specific values of the trans- 
port coefficients were then obtained by means of Fermi's 
golden rule, leading to a practically useful treatment. It 
is important to recognize that although the Langevin 
term may resemble the effect of the standard collision 
term, its origin is different: While the former results from 
the residual interaction, the latter arises from the quan- 
tum fluctuations that are inherent in the wave packets 
employed in the description. 

The proposed extension of the standard treatment rep- 
resents a formally well based approach towards incorpo- 
rating the effect of quantum fluctuations into the wave 
packet dynamics. Moreover, it leads to the desired sta- 
tistical properties in static scenarios and can be included 
in the dynamics in a conceptually simple and tractable 
manner. The method may therefore find useful appli- 
cation in the context of microscopic simulations of ac- 
tual many-body processes, such as fragment production 
in heavy-ion collisions. 
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APPENDIX A: MOMENTS OF THE 
TRANSITION RATE 

In this appendix, we show some technical details re- 
garding the calculation of the transition rate and its first 



and second moments which are required in sect. 44.3 



We start by rewriting the matrix element in the tran- 
sition rate J76|) , 

| < Z'\V\Z > | 2 = V(z', Z) V(Z, Z') I < Z'\Z > | 2 , 

(Al) 

where V{2',Z) = [Z'\V\Z)/(Z'\Z). The overlap be- 
tween two normalized states is approximately of Gaus- 
sian form, 



\<Z'\Z>\ 2 n exp (-SZ ■ C ■ SZ) , 
as can be seen by expanding in SZ = Z — Z, 
log | < Z'\Z > | 2 = -SZ -C-SZ + O ((Sz) 3 ) 



(A2) 



(A3) 



The above expression (A2) is consistent with the fact 
that dZ ■ C ■ dZ defines a infinitesimal squared dis- 
tance between two wave packets and, accordingly, that 
det C appears in the canonical measure. Furthermore, 
the relation is exact when C is unity, i.e. when anti- 
symmetrization is ignored. 

Using the expression (|A2|) and assuming that Z' E is 



close to Z E , we obtain 



< z'\v\z f 



> 



\v(z E ,z E )\ 2 e -(z' E -z E yc{z' E -z E ) 



(A4) 
(A5) 



Finally, when Z' E ss Z E then is small and the dis- 
tortion can be ignored, Z' w Z' E . Thus we arrive at the 
total transition rate given in eq. Jsi"|). 
We have used the resolution of unity, 



dV < Z\Z' > < Z'\Z > 



dV e -^-C-sZ 
I . (A6) 



and the first and second moment can be calculated by 
invoking the formulas 



J dV Sz n exp (SZ ■ C • SZ) = , 

dT' 5z n 8z n , exp (SZ-C ■ SZ) = (C -1 ) nn / 



(A7) 
(A8) 



In our previous work [|l5fl , we have subtracted the self- 
transition matrix element and expanded V in SZ. While 
this subtraction reduces the total transition rate, the re- 
lation between the drift and diffusion coefficients (|71 
remains valid. 
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